Floating phase in a dissipative Josephson junction array 
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We consider dissipative quantum phase transitions in Josephson junction arrays and show that 
the disordered phase in this extended system can be viewed as an unusual floating phase in which 
the states of local (0 + l)-dimensional elements (single Josephson junctions) can slide past each 
other despite arbitrary range spatial couplings among them. The unusual character of the metal- 
superconductor quantum critical point can be tested by measurements of the current voltage char- 
acteristic. This may be the simplest and most natural example of a floating phase. 

PACS numbers: 



It has recently been recognized [l| that phases and 
phase transitions in a given dimension can be embed- 
ded in a higher dimensional manifold. One realization 
of this is in a stack of two-dimensional (2D) layers of 
coupled XY spins with interlayer couplings along the 
layer- normals |l| . For a suitably chosen (and fairly com- 
plicated) set of competing interlayer couplings involving 
higher order gradients, it was shown in Ref. |lj that the 
3D-system can exhibit 2D-behavior: power-law ordered 
phases in each layer below a characteristic temperature, 
Tkt i and a high-temperature disordered phase with ex- 
ponentially decaying correlations. The transition sep- 
arating the two phases remains 2D Kosterlitz-Thouless 
type. Below Tkt, the 2D power-law phases in each layer 
'slide' past one another (a "floating" phase) Q. Strik- 
ingly, such a phase is impossible for more conventional 
Josephson-type couplings between the layers. Another 
interesting case is the recently discussed quantum theory 
of the smectic metal state in stripe phases , which also 
involves gradient couplings. Unfortunately, there are no 
definitive experimental realizations of a floating phase, 
although there is some evidence for it in lipid-DNA com- 
plexes |3| in which lipid bilayers are able to slide over 
each other without cost in energy . 

In this Letter we show that such a floating phase nat- 
urally occurs in the ground state of a system of coupled 
rcsistively shunted Josephson junctions [5|- The floating 
phase in this case is a state where the ground states of 
(0 + l)-dimensional elements (single Josephson junctions) 
can slide past each other despite coupling between them. 
Experimentally, this phase will behave like a metal with 
finite resistance^ Q- The striking point is that no addi- 
tional, higher order gradient couplings need be invoked, 
in contrast to the stack of XY spins treated in Ref. . 

The existence of such a floating phase in this system 
has been conjectured before (EH)!!- The arguments pre- 
sented there are, however, incomplete, since further than 
nearest-neighbor interactions, which are always present 
in any real system, are not considered in these treat- 



ments. This is particularly worrisome since the treat- 
ment of XY stacks found many regions of the phase 
diagram in which such interactions destroyed the float- 
ing phase. In this paper, we show that no such problem 
arises here: at least for sufficiently small further neighbor 
interactions, the floating phase is robust. 

A concrete realization of longer ranged Josephson cou- 
plings relevant to the present model is a system of super- 
conducting grains embedded in a metallic medium. The 
Josephson coupling Jy between the grains at rt and rj 
falls off very slowly at T = 0, specifically Jy ex 
in two dimensions |lfj||: this is due to pair field-pair field 
correlation function in a normal metal. 

The prior work also argued^, 0, 0| that there will be 
a phase transition as the shunt resistance is varied, be- 
tween the floating phase and a (Z?+l)-dimensional super- 
conducting phase, in which the phases of the junctions 
are coherently coupled. In addition to confirming this, 
we show from a renormalization group analysis that this 
transition, for suitable shunt resistance and small nearest 
neighbor Josephson coupling V, is governed by a line of 
quantum critical points with continuously varying crit- 
ical exponents, which in dimensions higher than one is 
very rare indeed. As we discuss below, the upper critical 
dimension in this problem is infinity. The phase diagram 
is shown in Fig. QJ In this figure, a = R/Rq, where R is 
the shunt resistance, and Rq = h/Ae 2 is the quantum of 
resistance. 

A similar phase diagram has been seen in a recent set of 
experiments ^]] , but the existence of continuously vary- 
ing critical exponents is not yet established. We em- 
phasize, however, that the floating phase is a metal, not 
an insulator: the metallic shunts short the current when 
the phase difference between the junctions become in- 
coherent. The metallic behavior was not established in 
this experiment, but only the insulating behavior of the 
unshunted array. Clearly, further experiments at suffi- 
ciently low temperatures will be necessary in establish- 
ing the true nature of the floating phase. We expect that 
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FIG. 1: Phase diagram and perturbative RG flows for ID, 
2D triangular lattice (3D face centered cubic), and 2D square 
lattice (3D simple cubic). M and SC indicate metal and su- 
perconductor, respectively. The thick lines correspond to line 
of fixed points. Crosses indicate special points on the a — 
axis denoting (D + 1) -dimensional XY transition. 



the current-voltage characteristic in the floating phase to 
be a non-universal power law controlled by dissipation 
a, similar to the results obtained previously |12(. In the 
superconducting state, the temperature-dependent non- 
universal power law will be similar to the classical XY 
model due to vortex unbinding in the presence of an im- 
posed current [l3| . 

As mentioned earlier, the phase transition is very 
strange, exhibiting continuously varying, non-universal 
critical exponents that depend on where the phase 
boundary is crossed. Additionally, the universality class 
of the transition depends not at all on the dimensional- 
ity of the system, but only on the local topology of the 
lattice. In particular, for a one-dimensional chain, the 
phase boundary in the V-a plane is vertical for small V 
(that is, independent of V), while for any lattice of any 
dimensionality that has closed loops of 3 nearest neigh- 
bor bonds (e.g., a 2D triangular or a 3D FCC lattice), 
the critical Josphson coupling V c obeys, for a near a c : 



V c (a) = A(a c - a), 



(1) 



where A is a non-universal constant. For any lattice of 
any dimensionality in which the smallest closed loops 
have more than 3 nearest neighbor bonds (e.g., a 2D 
square or 3D simple cubic lattice), the critical Josphson 
coupling V c obeys, for a near a c : 



V c (a) = y/A'{a c -a), 



(2) 



where A 1 is another non- universal constant. We will here- 
after refer to these distinct lattice types as "triangular- 
type" and "square-type" respectively. 

One non-universal critical exponent that can easily 
be measured experimentally is the thermal exponent v, 
which controls the scaling of characteristic temperatures 
T c h with distance from the T — quantum phase transi- 
tion. One such characteristic temperature for a's and V's 
that lie in the superconducting phase at T — is super- 
conducting to normal transition temperature T c . For a's 



and V's that lie in the floating phase at T = 0, a char- 
acteristic temperature is the temperature at which the 
resistivity of the array as a whole has its minimum Q. 
We predict that 



T ch cx (8V) V , 



(3) 



with SV = \V — V c \ . With a* ^ a c the value of a at 
which the transition occurs, the non-universal exponent 
v(a*) is given by 



v{a*) 



1 



A"(a c -a*Y 



(4) 



with A" = 1(2) for "triangular (square)-type" lattices. 

The array of the Josephson-junctions coupled to lo- 
cal Ohmic dissipation is described by the following 
imaginary-time (r) action: 

p@ i fin 

s/h = I ^w ^ { ~d )2dT+v ^ [1 ~ cosA9ij{T)]]d7 



(5) 



Here the sum is over nearest neighbor pairs in a 

lattice of arbitrary dimensionality and type. A8ij(u) n ) 
is a Fourier component of the phase difference be- 
tween the two superconducting grains, A^y(r), with 
uj n = ±27rn//3: n is an integer, and (3 is the inverse- 
temperature. The charging energy of a grain, Eq, is 
given by Eq cx 1/C, where C is the capacitance of the 
grain. The term containing a arises from integrating out 
the degrees of freedom of the harmonic-oscillator bath 
with a linear spectral-function (Ohmic dissipation) |l5j| . 
An Ohmic shunt, by definition, transfers charge contin- 
uously, while in quasiparticle tunneling in a supercon- 
ductor the charge is transferred in discrete quantized 
amounts [lr3 |. 

The term, Sj, describes inter-junction couplings medi- 
ated by longer-range Josephson-interactions among the 
grains, which has not been considered previously. These 
couplings are described by the general action, 



Sj/h. 



[Si] 



d,Tj([si\) cos [2J Sj9(fi 



■S S ,r)] 



(6) 

Here, 9[fi, r)'s are the phases of the order parameters in 
the superconducting grains, the S^s are arbitrary sep- 
aration vectors on the lattice (to allow for arbitrary 
ranged couplings) and Sj is an integer-valued function 
of the layer number i satisfying s; = 0. This last 
constraint follows from the absence of an external field, 
which implies "rotation invariance" under adding the 
same constant to all of the 6's. Note that the special 
case so = +l,^o = 0, s± = — l,Si = a 7 , where a 7 is 
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a nearest neighbor vector, is just the nearest neighbor 
Josephson coupling V in Eq. JSJ). 

To establish the existence of a dissipation-tuned quan- 
tum phase transition, we use a renormalization group 
analysis that is perturbative in V and J. We divide 
the field 9i(r) into slow and fast components 6i S (r) and 
Oi f (t) and write the partition function Z as 



Z = Z /J] m is{r) exp [-S s /h + ln(e 



-S'/h\ 



)o f (7) 



Here Zq is a normalization constant, Sq is the slow- 
frequency component of the quadratic part of the action 
containing the first and the third terms in Eq. (J5J, S' 
contains the remaining two terms, and (. • .)o/ denotes 
the average with respect to the fast components of the 
quadratic part. The average in Eq. Q contains a product 
of two factors involving V and J. 

To leading order in V and J, the coupling constants 
do not mix. After computing the averages, we rescale r, 
t' = r/b, b being the scale factor, to restore the original 
frequency cut-off, and then redefine the coupling con- 
stants to complete the renormalization. The dissipation 
term in Eq. JSJ is dimensionless, and so it is held fixed by 
the RG. The first term then has r-dimension — 1. Thus, 
it is irrelevant in the RG sense, and will henceforth be 
dropped for the remainder of this paper [l7j. 

Writing b — e , where I > is infinitesimal, and zq — 
z/2, where z is the coordination number of the lattice, 
we have to one loop order ||: 



f -('--)" 

dl zga 



(8) 



In the absence of the J's, the a-axis with V = is a line 
of stable fixed points for a < 1/zq and a line of unstable 
fixed points for a > 1/zq. The physical interpretation of 
the stable fixed points is as follows: for a < 1/zq, the 
barrier V among the different potential minima of the 
variable Af^ is irrelevant, and the system described by 
A9ij exhibits local quantum fluctuations. 

This local criticality can be seen from the equal-time 
correlation function of the phase differences: 



(exp[igA%(r)]) oc l/(/3fiw c )"W 



(9) 



with 77(g) =™|^ r l'( ( l~ n ) 2 j v' = X/(z n a), q a real number, 
and hw c = Eq. Note that the simplest such correlation 
function one could imagine, namely q = 1, does not de- 
cay algebraically because = 0. This result for 77(g) is 
familiar for surface roughening [Til Il9| . A similar calcu- 
lation shows that the unequal time correlation functions 
of the bond variables are also critical: 



(ex P [ l(? (A%(T)-A%(0))]) = 



1 



(cj c t) 2 "(9) 



(10) 



Although the correlation functions of the bond variables 
algebraically decay with time in any dimension, the same 



is not true for the correlation of the site variables. We 
have computed the unequal time site correlation func- 
tions [13] and have found that they vanish as r — > 00; the 
power law form is obtained for only D > 3. The locality 
of this phase follows from the result that [2(|, at zero 
temperature, spatial correlation functions of the grains, 
< exp[iq(6i(T) — 0i +r (r))] >, are identically zero in any 
dimension; spatially separated grains are completely un- 
corrected. Each junction is a dynamical system on its 
own, and there are only short-ranged correlations be- 
tween spatially separated junctions. 

For a > 1/zq, on the other hand, the barrier V grows, 
and the individual quantum systems saturate in one of 
the potential minima. V is like a local field term for 
the 'bond-spin' variable which becomes relevant above 
a = 1/zq, and orders the individual bond-spins. The 
transition at a = X/zq is purely local. 

What happens when longer-ranged spatial couplings 
among them are introduced? Calculations precisely anal- 
ogous to that leading to Eq. |(HJ| give, to first order in J, 



dJ 



r(N) 



with 



(11) 



(12) 



1,3 



The "potential" U(r) = % £-(e^ - with fa = 

53 7 (1 — e lq ' a ~<), where the sum on 7 is over all nearest 
neighbors. It is straightforward to show that U(r) is 
nothing but the "lattice Coulomb potential" of a unit 
negative charge at the origin, with the zero of the po- 
tential set at r — 0. That is, U(r) satisfies the "lattice 
Poisson equation" : ^ 7 U(r — a 7 ) — zU(r) = —8?$- For a 
symmetrical (e.g., square, hexagonal, cubic) lattice, the 
left-hand side is just the "lattice Laplacian" , approaching 
V 2 U x 0(o 2 ) where (a = |a 7 |). 

The quantity T in Ullfl and H12|) is then clearly just 
equal to the potential energy of a neutral (since J2i s i — 
0) plasma of quantized (since the s^s arc integers) 
charges Si on the lattice. The most relevant J([sj]) 
is therefore clearly the one that corresponds, in this 
Coulomb analogy, to the lowest interaction energy. (Note 
that strictly speaking T corresponds to twice this energy, 
because the sum in l)12|l double counts). Aside from the 
trivial configuration in which all the Si = , the low- 
est energy configuration is clearly one in which there are 
two equal and opposite unit magnitude charges on near- 
est neighbor sites: i.e., sq = +1, So = 0, s% = —1, Si = a~ ( . 
As discussed earlier, this corresponds to the nearest- 
neighbor Josephson coupling in equation (0. Thus, we 
have established that that coupling is, indeed, the most 
relevant, as asserted earlier. Furthermore, using simple 
symmetry arguments, one can show that for a symmetric 
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lattice (e.g., square, hexagonal, cubic), where all nearest- 
neighbor sites are equivalent, t/(a 7 ) = — which re- 
covers the recursion relation JSJ for V. Hence, all other 
couplings are irrelevant for a < a c = l/z ; as a result 
they affect neither the floating phase nor the transition 
between it and the (D + l)-dimensionally coupled phase. 

In summary, as long as a < 1/zq, V, and all of the 
other J([sj])'s are irrelevant. The (0+l)-dimensional sys- 
tems constitute locally critical power-law phases. Since 
J's are also irrelevant in this regime, the local systems 
are in a floating phase [ij ■ As soon as V becomes relevant 
at a > 1/zq, it grows to oo. In this limit, this enormous 
periodic potential obviously traps the field A6 at the bot- 
tom of one of the minima of the periodic potential; and 
provides a mass for fluctuations about that minimum. 
Such a mass destroys the logarithmically divergent fluc- 
tuations, which give rise to the r([sj])/a term in the 
recursion relation in Eq. ^2 f° r J- Thus, in this limit 
that recursion relation just becomes dJ/dl = J. This 
in turn makes all J's relevant simultaneously with the 
nearest neighbor coupling V, ensuring that there are no 
further transitions as a is increased. Thus, there is only 
one global transition with a as the tuning-parameter. 

Even though the J's become relevant as soon as V 
does, they remain irrelevant far enough on the a-axis 
for the (0 + l)-dimensional transition to happen. The 
situation is reversed for planes of power-law XF-phases 
with interlayer couplings, and so Ref. had to intro- 
duce extra gradient-couplings among the planes to get 
floating power-law phases. In the present case, the metal 
to superconductor quantum phase transition obtained by 
tuning dissipation in arbitrary spatial dimension is then a 
true (0 + l)-dimensional phase transition. In the floating 
local critical phase on the metallic side of the transition, 
the local systems are, however, not decoupled; J's are ir- 
relevant, but not zero, and affect the physical properties 
of the metallic phase only perturbatively. 

To obtain higher order corrections to the recursion re- 
lation for V, Eq. |SJ we perform a cumulant expansion of 
Eq. in second order in V. We find no contribution to 
Eq. [HI in D — 1, while for D = 2 triangular lattice there 
is a contribution at order V 2 , which is 

) + c!v\ as) 

al zoa 

where C\ is a positive constant |2J|. This holds for any 
lattice in any dimension, as long as there are just three 
sides in a minimum closed loop. For 2D square lattice, 
however, there is no renormalization of V at order V 2 . 
Indeed, the first correction to Eq. [S] comes only at order 
V 3 |20| for any lattice where the minimum closed loop 
has number of sides more than three, 

^ = ni-— )+c 2 v 3 , (w) 

al zoa 
where C2 is another positive constant |20| . It is impor- 



tant to emphasize that the line of fixed points implied by 
Eq. US and [H exist in any dimension, which implies that 
the upper critical dimension is infinity. Combined with 
the flow equation for a, ^ — 0, which remains true in 
our perturbative treatment at all orders in V, we get the 
RG flows, the phase diagram, and the exponents quoted 
earlier. 

Not only the existence of a floating phase and its val- 
idation in the resistively shunted Josephson junction ar- 
rays are exciting possibilities, but the present analysis 
should find its use in a much broader context: local quan- 
tum critical points are widely discussed in heavy fermion 
systems [2lJ and are also argued to be crucial to the ubiq- 
uitous metallic phases observed in a diverse class of con- 
densed matter systems |22j |. 
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